1 



Second-order Gauge Invariant Perturbation Theory 

Perturbative curvatures in the two-parameter case 

Kouji Nakamura 

Department of Astronomical Science, the Graduate University for Advanced 
\Q • Studies, Mitaka, Tokyo 181-8588, Japan. 

O ' 

O ■ (Received October 6, 2004) 

(N 

Based on the gauge invariant variables proposed in our previous paper [K. Nakamura, 
Prog. Theor. Phys. 110 (2003), 723], some formulae for the perturbative curvatures of each 
order are derived. We follow the general framework of the second-order gauge invariant 
perturbation theory on an arbitrary background spacetime to derive these formulae. It is 
found that these perturbative curvatures have the same form as those given in the definitions 
of gauge invariant variables for arbitrary perturbatives fields, which were proposed in the 
above paper. As a result, we explicitly see that any perturbative Einstein equation can be 
given in terms of a gauge invariant form. We briefly discuss physical situations to which this 
framework should be applied. 
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In many theories of physics, realistic situations are often difficult to describe by 
an exact solution of a theory, because theories in physics and their exact solutions are 
often too idealized to properly represent natural phenomena. Given this situation, we 



q-( have to consider perturbative approaches to investigate realistic situations. General 

relativity is one theory in which the construction of exact solutions is not so easy. 
Though there are many exact solutions to the Einstein equation, 1 ) these are often 
too idealized. For this reason, general relativistic perturbation theory is a useful 
technique to investigate natural phenomena. 2 )' 3 ^ 

In addition to this technical problem, general relativity is based on the concept 
of general covariance. Intuitively, the principle of general covariance states that 
there is no preferred coordinate system in nature, though the notion of general 
covariance is mathematically included in the definition of a spacetime manifold in a 
trivial way. This is based on the philosophy that coordinate systems are originally 
chosen by us, and natural phenomena have nothing to do with our coordinate system. 
Due to this general covariance, the gauge degree of freedom, which is an unphysical 
degree of freedom of perturbations, arises in general relativistic perturbations. To 
obtain physically meaningful results, we have to fix this gauge degrees of freedom 
or to extract the gauge invariant part of perturbations. A similar situation has been 
found in recent investigations of the oscillatory behavior of a gravitating Nambu- 
Goto membrane, 4 ^' 5 ) which concern the dynamical degrees of freedom of extended 
gravitating objects. 

On the other hand, higher-order multi-parameter perturbations, in which there 
are two or more small parameters, can be applied to many physical situations. One 
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well-known application of two-parameter perturbation theory is perturbations of 
a spherical star, 6 ) in which we choose the gravitational field of the spherical star 
as the background spacetime for the perturbations, one of the parameters for the 
perturbations corresponds to the rotation of the star, and the other is its pulsation 
amplitude. The effects due to the rotation-pulsation coupling are described at higher 
orders. A similar perturbation theory on the Minkowski background spacetime was 
developed by the present author to study the comparison of the oscillatory behav- 
ior of a gravitating string with that of a test string. 5 ) Even in the one-parameter 
case, it is interesting to consider higher-order perturbations. In particular, Gleiser 
et al. 7 ) reported that second-order perturbations yield accurate wave forms of grav- 
itational waves. Hence, it is worthwhile to investigate higher-order multi-parameter 
perturbation theory from a general point of view. 

Motivated by these physical applications, the general relativistic gauge invari- 
ant multi-parameter perturbation theory has been developed in a number of pa- 
pers. 5 )' 8 )' 9 ) In particular, the procedure to find gauge invariant variables for higher- 
order perturbations on a generic background spacetime was proposed by the present 
author 8 ^ assuming that we have already known the procedure to find gauge invari- 
ant variables for a linear-order metric perturbations. The contents of this paper are 
based on this proposal. The main purpose of this paper is to present some formulae 
of the second-order perturbative curvatures within the two-parameter perturbation 
theory that are useful in some physical applications. When we derive these formulae, 
we follow the general framework of the second-order gauge invariant perturbation 
theory on an arbitrary background spacetime. This framework is originally proposed 
by Stewart et al. 10 ) and developed by Bruni et al., 9 )' 11 ) and the present author. 8 ) 
These perturbative curvatures have the same form as those given in the definitions 
of gauge invariant variables for arbitrary perturbative fields which are proposed in 
a previous paper. 8 ) As in that paper, we do not make any specific assumption re- 
garding the background spacetime and the physical meaning of the two-parameter 
family. Because we make no assumption concerning the background spacetime, this 
framework has a wide area of applications. 

The organization of this paper is as follows. In we review the general frame- 
work of the second order gauge invariant perturbation theory. We mainly review 
the one-parameter perturbation theory. We emphasize that the review in 32 of this 
paper is based on the idea of Stewart et al. 10 ) and the development carried out by 
Bruni et al. 9 )' 11 ) In we present some formulae for the second-order perturbative 
curvatures within the two-parameter perturbation theory. We also give a derivation 
of these formulae. The final section, is devoted to a summary and brief discussion 
of physical situations to which this framework of higher-order perturbation theory 
should be applied. We employ the notation of our previous paper 8 ) and use the 
abstract index notation. 12 ) 

§2. Gauge degree of freedom in perturbation theory 

In this section, we briefly review the gauge degree of freedom in general rela- 
tivistic perturbations. This was originally discussed by Stewart et al. 10 ) To explain 
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the gauge degree of freedom in perturbation theories, we have to keep in mind what 
we are doing when we consider perturbations. We first comment on the intuitive ex- 
planation of the gauge degree of freedom in H2.ll Next, in N2.21 we review the more 
precise mathematical formulation of the perturbations in the theories with general 
covariance. The explanation given here is based on the works of Bruni et al., 11 ) which 
represent extensions of the idea of Stewart et al.. When we consider the perturba- 
tions in the theory with general covariance, we have to exclude these gauge degrees 
of freedom in perturbations. To accomplish this, gauge invariant quantities of the 
perturbations are useful, and these are regarded as physically meaningful quantities. 
In ^2.31 based on the mathematical preparation given in ^12.21 we review the proce- 
dure to find gauge invariant quantities of perturbations, which was developed by the 
present author. 8 ) 

2.1. What is "gauge"? 

In perturbation theory, we always 
treat two spacetime manifolds. One is 
the physical spacetime, which we at- 
tempt to describe by perturbations, and 
the other is the background spacetime, 
which we prepare for perturbative anal- 
yses. Let us denote the physical space- 
time by (Ai,g a b) an d the background 
spacetime by (Mo,g a b)- Keeping in 
mind these two spacetime manifolds, let 
us formally denote the spacetime metric 
and the other physical tensor fields on 
the physical spacetime M by Q. As the perturbation of the physical variable Q, we 
always write equations of the form 




Fig. 1. "gauge choice" in one-parameter per- 
turbation. 



Q( i y) = Qo(p) + SQ(p). 



(2-1) 



Usually, this equation is regarded as the relation between the physical variable Q 
and its background value Qq of the same field, or simply as the definition of the 
deviation 5Q of Q from its background value, Qq. In fact, through the equation 
(|2-1[) . we have implicitly assigned a correspondence between points of the physical 
and the background spacetime, since this equation gives a relation between field 
variables Q, Qo and 6Q. More specifically, Q("p") in the left-hand side of Eq. 1)2-1(1 
is a field on the physical spacetime M, and "p" E M. On the other hand, we should 
regard the background value Qo(p) of Q( a p") and its deviation 5Q(p) from Qo(p), 
which are on the right-hand side of Eq. (|2-1|) , as fields on the background spacetime 
Mo and p € Mo- Because Eq. (|2-ip is for field variables, it implicitly states that 
the points "p" £ M and p 6 Mq are same. This is an implicit assumption of the 
existence of a map A^o M : p € Mo i— > "p" € M, which is usually called a "gauge 
choice" in perturbation theory. 10 ) Clearly, this is more than the usual assignment of 
coordinate labels to points within the single spacetime. 

It is important to note that the correspondence between points on each M e , 
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which established by such a relation as Eq. (|2-ljl . is not unique to the theory in 
which general covariance is imposed. Rather, Eq. (|2- 1|) involves the degree of freedom 
corresponding to the choice of the map X : Mq —* A4 (the choice of the point 
identification map — > «M). This is called the gauge degree of freedom. Further, 
such a degree of freedom always exists in the perturbations of a theory in which we 
impose general covariance, unless, there is a preferred coordinate system in the theory 
and we naturally introduce this coordinate system onto both Mo and M.. Then, we 
can choose the identification map X using this coordinate system. However, there 
is no such coordinate system, due to general covariance, and we have no guiding 
principle to choose the identification map X. Therefore, we may identify "p" G M 
with q G M.Q (q ^ p) instead of p G M.$. (See Fig. ^) A gauge transformation is 
simply a change of the map X. 

2.2. More precise formulation of perturbations 

In this section, we review the more precise formulation concerning about " 'gauge 
degree of freedom! 7 based on the above understanding of "gauges" , 10 )' 11 ) We mainly 
review the one-parameter perturbation theory in H2. 2. H and comment on the results 
in two-parameter perturbation theory in ^2,2.21 The essential part of the multi- 
parameter perturbations is completely similar to the one-parameter case. 8 )' 11 ) De- 
tails can be seen in the papers by Bruni et al. 11 ) and by the present author. 8 ) 

2.2.1. One-parameter perturbation theory 

We denote the perturbation parameter by e, and we consider the m+l-dimensional 
manifold J\f = Ai x R, where m = dim M and e G M, as depicted in Fig. ^ By 
this construction, the manifold J\f is foliated by m-dimensional submanifolds Ai t 
that are diffeomorphic to the physical spacetime M.. The background M.q = A/"| e=0 
and the physical spacetime M. = A4 t = AA| R=e are also submanifolds embedded in 
N. Each point on M is identified by a pair [p, e), where p G A4 e , and each point 
on the background spacetime .Mo m N is identified by e = 0. The manifold M 
has a natural differentiable structure consisting of the direct product of M and K. 
By this construction, the perturbed spacetimes A4 e for each e must have the same 
differential structure. In other words, we require that perturbations be continuous 
in the sense that (M.,g a b) and (M.o,g a b) are connected by a continuous curve within 
the extended spacetime N '. Hence, the changes of the differential structure resulting 
from the perturbation, for example the formation of singularities and singular per- 
turbations in the sense of fluid mechanics, are excluded from the study carried out 
in this paper. 

Let us consider the set of field equations 

£[Qe] = (2-2) 

on the physical spacetime M. e for the physical variables Q t on A4 e . The field equa- 
tion H2-2|) formally represents the Einstein equation for the metric on A4 e and the 
equations for matter fields on Ai e - If a tensor field Q e is given on each Ai £ , Q e is 
automatically extended to a tensor field on N by Q(p, e) := Q e (p), where p G M. e . In 
this extension, the field equation (|2-2|) is regarded as the equation on the extended 
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manifold TV. Thus, we have extended an arbitrary tensor field and field equations 
(|2-2[) on each A4 e to those on the extended manifold Af. 

Tensor fields on Af obtained by the above construction are necessarily "tangent" 
to each A4 t , i.e., their normal component to each A4 e identically vanishes. To con- 
sider the basis of the tangent space of Af, we introduce the normal form of each A4 e 
in Af and its dual. These are denoted by (de) a and (d/de) a , respectively, and they 
satisfy 

d 



(& H*J =L (2 ' 3) 
The form (de) a and its dual (d/de) a are normal to any tensor field that is extended 
from the tangent space on each A4 e by the above construction. The set consisting 
of (de) a , (d/de) a and the basis of the tangent space on each AA e is regarded as the 
basis of the tangent space of Af. 

To define the perturbation of an arbitrary tensor field Q, we compare Q on the 
physical spacetime A4 e with Qq on the background spacetime, and it is necessary 
to identify the points of A4 t with those of A4q. This point identification map is the 
so-called gauge choice in the context of perturbation theories, as mentioned in ^2.11 
The gauge choice is accomplished by assigning a diffeomorphism X e : Af — > Af such 
that X e : A4.Q — > A4 e . Following the paper of Bruni et al., 11 ) we introduce a gauge 
choice X t as one of the one-parameter groups of diffeomorphisms that satisfy the 
property 

X ei +e 2 = X 62 = X t2 o X ei . (2"4) 
This one-parameter group of diffeomorphisms is generated by the vector field X f)? e y 
This vector field X ffi e \, which we call a generator, is defined by the action of the 
corresponding pull-back X* for a generic tensor field Q on A4 x M : 



X*Q - Q 



£^ (e) Q:=lim ^ (2-5) 



and it is decomposed as 



%■) =: + a , 9 a (de) a = 0, £e_0 a = 0. (2-6) 

The third condition in (|2-6(l is imposed merely for simplicity. Except for the condi- 
tions in ()2-6|) . we may regard 6 a as an arbitrary vector field on A4 t (not on Af); i.e., 
the arbitrariness of the gauge choice is given by that of the vector field 9 a . 
The Taylor expansion of the pull-back X*Q is given by 



00 f k T r) k 1 °° f k 

fc=0 L J e=0 fc=0 



Once the definition of the pull-back of the gauge choice X t is given, the perturbation 
A x Q e of a tensor field Q under the gauge choice X e is simply defined as 

A x Q e := X:Q\ Mq - Q . (2-8) 
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We note that all the variables in this definition are defined on .Mo- The first term 
on the right-hand side of (|2-8j) can be Taylor-expanded as 

OO f, 

k=0 



Equations (|2-8|) and (|2-9|) define the perturbation of 0(k) of a physical variable Q e 
under the gauge choice X and its background value #Q = Qq. Through Eqs. H2-7j) 
and (|2-9|l . each order perturbation X Q under the gauge choice X e is given by 

( 5Q = £\Q KA ■ (2-io) 

Mo 

The above understanding of the gauge choice and perturbations naturally leads 
to the gauge transformation rules between different gauge choices and the concept 
of gauge invariance as follows. 

Suppose that X e and y e are two one-parameter groups of diffeomorphisms with 
the generators x rf and ^if on jV, respectively, i.e., X e and y e are two gauge choices. 
These generators are decomposed in the same manner as Eqs. (j2-6j) : 

x v a = (^J + o a , y v a = ^ a + c a . (2-n) 

The integral curves of each x n a and ^rf 1 in M are the orbit of the action of the gauge 
choices X e and y e , respectively. Since the generators x rf and ^77° are transverse to 
each M e everywhere on J\[, the integral curves of these vector field intersect with 
each Ai E - Therefore, points lying on the same integral curve of either of the two are 
to be regarded as the same point within the respective gauges (see Fig.^). Hence, X t 
and y t are both point identification maps. When 9 a 7^ i a , these point identification 
maps are regarded as two different gauge choices. 

Suppose that X e and y e are two different gauge choices which are generated by 
the vector fields x rf and ^rf, respectively. These gauge choices also pull back a 
generic tensor field Q on N to two other tensor fields, X*Q and y*Q, for any given 
value of e. In particular, on .Mo, we now have three tensor fields associated with a 
tensor field Q; i.e., one is the background value Qq of Q, and the other two are the 
pulled back variables of Q from A4 e to Mq by the two different gauge choices, 



00 k 

Q e := x:q\ Mq = ~n xQ = Qo + A*Qe, (2-12) 



A;! 

k=0 



Qe ■■= y:Q\ Mo = E uy Q = Qo + aVq *- (2-13) 



k=0 



Here, we have used Eqs. (|2-8(l and (|2-9|) . Because X t and y t are gauge choices that 
map the background spacetime Mo into the physical spacetime M. e , x Q e and ^Q e 
are the representations on A4q of the perturbed tensor field Q in the two different 
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gauges. The quantities and ^yQ in Eqs. (|2-12j) and (|2T;-i|) are the perturbations 
of 0{k) in the gauges X and y, respectively. 

Now, we consider the concept of gauge invariance. Following the paper of Bruni 
et al., 9 ) we consider the concept of gauge invariance up to order n. We say that Q 
is gauge invariant up to order n iff for any two gauges X and y 

Wq = yk, with k < n. (2-14) 

From this definition, we can prove that the nth-order perturbation of a tensor field 

(k) 

Q is gauge invariant up to order n iff in a given gauge X we have £^ X Q = for any 
vector field £ a defined on Mo and for any k < n. As a consequence, the nth-order 
perturbation of a tensor field Q is gauge invariant up to order n iff Qo and all its 
perturbations of lower than nth order are, in any gauge, either vanishing or constant 
scalars, or a combination of Kronecker deltas with constant coefficients. 9 )~ n ) 

In general, the representation X Q € on Mo of the perturbed variable Q on M e 
depends on the gauge choice X t . If we apply a different gauge choice, the repre- 
sentation of Q e on Mo may change. Recalling that the gauge choice X is a point 
identification map from Mo to M e (see Fig.^), the change of the gauge choice from 
X e to y e is represented by the diffeomorphism 

$ e :=(X e )- 1 oy e . (2-15) 

This diffeomorphism <P e is the map <P t : Mo — * Mo for each value of e € R. As shown 
in Fig.^ the diffeomorphism <£ e changes the point identification, as expected from the 
understanding of the gauge choice discussed in ^12.11 Therefore, the diffeomorphism 
<P e is regarded as the gauge transformation ^> t : X t — > y e . 

The gauge transformation <? e induces a pull-back from the representation x Q e 
of the perturbed tensor field in the gauge choice X e to the representation ^Q e in the 
gauge choice 34- Actually, the tensor fields x Q e and y Q e , which are defined on Mo, 
are connected by the linear map <&* as 



y Qe = ytQ\ Mo = (ytix^YQ 



Mo 



= (X-^y e y(X:Q) = $**Q e . (2-16) 

Mo 

According to generic arguments concerning the Taylor expansion of the pull-back of 
a tensor field on the same manifold, 8 )' n ) the gauge transformation <P* x Q e should be 
given by the form 

K X Q = X Q + e£& X Q + J {£& + £&*Q + G(e 3 ), (247) 

where the vector fields £f and £f are the generators of the gauge transformation <P e . 

Comparing the representation (|2-17j) of the expansion in terms of the generators 
£p of the pull-back <P* X Q and that in terms of the generators X, rf? e \ and ^fff e \ of 

the pull-back y* o (A^ -1 )* Q (= <P* X Q), we easily find explicit expressions for the 
generators of the gauge transformation <P = X^ 1 o y in terms of the generators 
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V? € ) ana - ^(e) °^ ^ ne S au S e choices- Further, because the gauge transformation <P e 
is a map within the background spacetime Mo, the generator should be given as 
vector fields on A4q. The explicit expression of the generators £p in terms of the 
components of the generators of the gauge choices is given in some papers. 8 )' 11 - 1 

We can now derive the relation between the perturbations in the two different 
gauges. Up to second order, these relations are derived by substituting ()2T2j) and 
(t2~T31l into (I2~T7I) : 

yQ - ^xQ = ^(iQoi ( 2 '18) 
( y ] Q - ®Q = 2£ kl $Q + {£ k2) + Qo- (2-19) 

These results are, of course, consistent with the concept of gauge invariance up to 
order n, as introduced above. Inspecting these gauge transformation rules, we can 
define the gauge invariant variables. 

2.2.2. Two-parameter perturbation theory 

Here, we briefly review the two-parameter case. We denote the two parameters 
for the perturbation by e and A. In this case, we have to consider the extended 
manifold J\f = A4 x M 2 , instead of TV = M x R used in the one-parameter case, 
where (e,A) £ M 2 . As in the one-parameter case, the gauge choice X e ^\ is a point 
identification map X e \ : M.q —> M on M. This gauge choice X e \ has the property 

^ei.Ai ° ^2,A 2 = X ei+e2j \ 1+ \ 2 Ai, A2, ei, 62 € M. (2-20) 
This property implies that 

%e,\ = Xe,0 ° <^0,A = %0,\ ° <*e,0, (2-21) 

where X e> o and Xq^\ are two one-parameter groups of diffeomorphisms defined by the 
property Eq. (|2-4|) . 

We denote the generators of X e ^ and Xo t \ by ?7? e ) an d ^(A)' respectively. We 
also introduce the basis (d/de) a , (de) a , (d/d\) a , (d\) a as vector fields on M, which 
satisfy conditions similar to those in (|2-3|) . 8 ) Using these basis, the generators x rj? e \ 

( X V(\)) °f the one-parameter group of diffeomorphisms X e ^ (<^o,a) is decomposed in 
the same manner as in Eq. ()2-6|) . The property ()2-21|) is expressed by 

[%),^(A)] a = (2-22) 

in terms of these generators. The Taylor expansion of the pull-back X* X Q is given 
by 



°° \fc Jc 



klk'l 

k,k'=0 



ak+k' 
dX kQ e k> ^,A^ 



(2-23) 

A=e=0 



00 \kk' 
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The perturbation AqQ € \ of an arbitrary tensor field Q in terms of the gauge choice 
X e> \ is given by 

°° \ k f k ' (V v\ 

A$Q*,x ■■= <xQ\ Mo - Qo, K>Q\ Mo =Y.Tw ( 'xQ, ( 2 - 25 ) 

k,k'=0 



where ' -JQ is the perturbation of order (A;, fc') and ' = Qq. Together with the 

,k>) 
x 



expansion given in (|2-24l) and (|2-25|) . each order perturbation with the gauge 



choice X e> \ is given by 

(2-26) 



(k,k') n _ p k ok' n 



Mo 



The concept of different gauges, a gauge transformation, gauge invariance, and 
the definition of gauge transformation rules in the two-parameter case are similar to 
those in the one-parameter case. For this reason, we can derive the following gauge 
transformation rules: 8 ) 

{p iQ- (p iQ = £^ q) Qo for (p,g) = (l,0),(0,l), (2-27) 

{P 'yQ ~ ^' q X Q = 2£ £, { v^ 2 "xQ + I^Cp.?) + °^(§,§)} 

for (p,g) = (2,0),(0,2), (2-28) 

+ |^d,i) + 2^(1.0)^(0,1) + 2^(0,1)^(1,0)} Qo, (2-29) 

where the are the generators for gauge transformation <P €j \ := (X ty \) 1 o 34,A- 
In this paper, we treat these gauge transformation rules of two-parameter per- 
turbation theory as mentioned in the introduction (JQ), because the one-parameter 
case considered above can be treated as a special case of this two-parameter case. 

2.3. Gauge invariant variables 

Inspecting the gauge transformation rules (|2-27|) - (|2-29|) . we can define the gauge 
invariant variables for a metric perturbation and for arbitrary matter fields. Em- 
ploying the idea of gauge invariance up to order n for nth-order perturbations, n ) 
we proposed the procedure to construct gauge invariant variables of higher-order 
perturbations. 8 ^ This proposal is as follows. First, we construct gauge invariant 
variables for the metric perturbation. Then, we define the gauge invariant variables 
for perturbations of an arbitrary field, excluding perturbations of the metric. The 
procedure to find the gauge invariant part of a higher-order perturbation is a simple 
extension of that for linear-order perturbations. 

To consider the metric perturbation, we expand the metric on the physical space- 
time A4, which is pulled back to the background spacetime A4q using a gauge choice 
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in the form given in 1)2-2 



= £ wr (fe,fc X» (2-30) 

fc',fc=0 

. (1,0), . \ (0,1) T 

= 9ab + e ^aft + A' x h ah 

+^ {2 ' 0) X h ab + eX^hat + ^ m x h ah + O s (e, A), (2-31) 

where (°' )/i a ft = g a b is the metric on the background spacetime M®. Of course, the 
expansion 1)2-31)1 of the metric depends entierly on the gauge choice X. Nevertheless, 
we do not explicitly express the index of the gauge choice X in an expression if there 
is no possibility of confusion. 

Our starting point to construct gauge invariant variables is the assumption that 
we already know the procedure to find gauge invariant variables for the linear metric 
perturbations. Then, linear metric perturbations ^'°^h a b ft ' 1 ' hub) are decomposed 
as 

Mhab =: ^H ab + 2V {a Mx b) , (p,q) = (1,0), (0,1), (2-32) 

where ^H ab and Mx a are the gauge invariant and variant parts of the linear- 
order metric perturbations. 8 ) Hence, under the gauge transformation 1)2-27)1 . these 
are transformed as (p '$W afe - {p ' q) x H ab = and {p ' q jx a - M x X a = & y 

As emphasized in a previous paper, 8 ) the above assumption is quite strong and 
it is not trivial to carry out the systematic decomposition 1)2-32)1 on an arbitrary 
background spacetime, as this procedure depends completely on the background 
spacetime (Mo,g a b)- However, this procedure is known in the perturbation theory 
on some simple background spacetimes, for example the cosmological perturbations 
of homogeneous and isotropic universes 2 ) or perturbations of spherically symmetric 
spacetimes. 3 ^ Further, from a general point of view, knowledge of linear perturbation 
theory is always necessary to carry out the second-order perturbations. For these 
reasons, we start from this assumption in spite of the fact that it is quite strong. 

Once we accept this assumption, we can always find gauge invariant variables 
for higher-order perturbations. 8 ) As shown in a previous paper, 8 ) at second-order, 
the metric perturbations are decomposed as 

^h a b=- {j> ' q) n ah + 2£ m)x ^U) hah 

+ (£ Mx - £ 2 ( P 4)x )g ab , (p,g) = (2,0),(0,2); (2-33) 

^'^hab =: ^'^Hab + £ (.0,1) X^ ^ab + £ '(i,o) x^ ^ab 

+ i)x ~ 2^<- 1 ' 0) x£(.°Vx ~ 2"^ - 1 )X'Ai.o)x| 9ab, (2-34) 

where ^H ab and ^X a are the gauge invariant and variant parts of the metric 
perturbations under the gauge transformation rules (|2-27j) - l|2-29j) . 
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Furthermore, using the gauge variant parts (P> q >X a of metric perturbations, 8 ) 
gauge invariant variables for an arbitrary field Q excluding the metric are given by 8 ) 

MQ := - £ (p , q)x Q , (p, q) = (1, 0), (0, 1), (2-35) 

-{£ {P , q)x - £ 2 { p^ )x }Qo, (p,g) = (2,0),(0,2), (2-36) 

(i,D Q := (l,i) Q _ £(1 , 0)X (^)Q - £ (0 , 1)X (^Q 

- i)x ~ 2^( 1 '°)x^'( ' 1 )x ~~ 2'^ (0 ' 1) ^^' (1 ' 0) ^| < ^°" (2-37) 

It is straightforward to confirm that the variables defined by (|2-35|) - ()2-37j) are 

gauge invariant under the gauge transformation rules (j2-27j) - l|2-29|) . In this paper, 
we derive some formulae for second-order perturbations of curvatures from the ex- 
pansion of the metric perturbation on a generic spacetime. The starting point of this 
derivation is the decomposition (|2-32|) - (|2-34j) of the metric perturbations in terms of 
the gauge invariant and variant variables. As a result, we find that all formulae have 
forms which are similar to those given in the definitions (|2-35|) - (|2-37|) of the gauge 
invariant variables for arbitrary matter fields. 



§3. Formulae of perturbative curvatures 



Now, we derive the formulae for the perturbative curvatures at each order in two 
parameter perturbation theory, following the standard derivation of the perturbative 
curvature. 12 ^ 

The starting point of the derivation is simply the definition of the curvature 
R abc d on the physical spacetime (M,g ab ) 

(V„ V b - V 6 V ) uo c = u> d R abc d , (3-1) 

where V a is the covariant derivative associated with the metric g ao on the physical 
spacetime M and u> c is an arbitrary one- form on the physical spacetime M. We 
similarly define the curvature R abc d on the background spacetime {M.o,g a b)i 

(V a Vb ~ V 6 V a ) uj c = uj d R abc d , (3-2) 

where V a is the covariant derivative associated with the metric g ab on the background 
spacetime Mo, and ui c is an arbitrary one- form on the background spacetime Mo. 
Our task is to compare R abc d and R abc d - To accomplish this, we have to consider 
the difference between the gauge choices for the physical spacetime M and the 
background spacetime Mq as discussed above. 

To compare the Riemann curvature (|3T[) of the physical spacetime M and 
that (|3-2|) of the background spacetime Mo, we introduce the derivative operator 
X* V a on the background spacetime A^o- This derivative operator X* V a (X^ 1 ) 

is the pull-back of the covariant derivative V a associated with the metric g ab on the 
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physical spacetime M.. The property of the derivative operator X*X7 a (X x )* as the 
covariant derivative is given by 

X*V a ({X- 1 )* X*9ab) =0, (3-3) 

where X*g ab is the pull-back of the metric on the physical spacetime M, which is 
expanded as Eq. (|2-31|) . Through the introduction of this operator X*V ' a (<-f _1 ) , we 
can regard the definition of the Riemann curvature (|3-1|) on the physical spacetime 
M. as an equation on the background spacetime. Since the pull-back X*V a (X^ 1 ) 
on the background spacetime A4q of the covariant derivative V a on the physical 
spacetime M is linear, satisfies the Leibnitz rule, commutes with contraction, is 
consistent with the concept of tangent vectors, and is torsion free, *) we can regard 
it as a derivative operator on the background spacetime .Mo- 12 '* Of course, the 
representation of this derivative operator X*V a (X" 1 ) on the background spacetime 
M.0 depends entirely on the gauge choice X. Though we should keep in mind that 
we have already chosen a gauge when we regard Eq. (|3Tj) as an equation on the 
background spacetime M.q, we do not explicitly express the index of the gauge choice 
X in any expression, again. 

Since V a (= X*V a (<^ -1 ) ) may be regarded as a derivative operator on the 
background spacetime that satisfies V a g bc = 0, there exists a tensor field C c ab on 
the background spacetime Mq such that 

V a w 6 = V a oj b - C c ab uj c , (3-4) 

where u) c is an arbitrary one- form on the background spacetime Mq. From the 
property (|3-3j) of the covariant derivative operator V a on M., the tensor field C c ab 
on .Mo is given by 

C c ab = \g cd (V a g db + V b g da - V d g ab ) . (3-5) 

We note that the gauge dependence of the derivative V a as an operator on A4$ 
is included only in this tensor field C c ab . From Eq. (|3-1|) . the Riemann curvature 
R abc d associated with the metric g ab is given by the Riemann curvature R abc d on 
the background spacetime and the tensor field C c ab as follows: 

Rabc d = R-abc d ~ ^[a^ d b}c + ^'cla^V' 

To obtain a perturbative expression of the curvatures, we first calculate the expansion 
of the inverse metric g ab , and then the perturbative expression of the tensor C c ab by 
using Eq. (|3-5(l . Next, we derive an expression of the perturbative curvature. 

In this paper, we present some formulae for the second-order perturbative curva- 
ture in the two-parameter perturbation theory. To derive the second-order formulae, 
we first calculate the 0(e\) formulae, since the other second-order formulae [0(e 2 ) 
and 0(A 2 )] are easily derived from these for O(eA) through a simple replacement of 

*' In this paper, we do not treat the torsion tensor. If we wish to consider a spacetime with 
torsion, we have to extend the formulation to that including the torsion tensor. 
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the perturbative variables. We also note that all variables on the physical spacetime 
M. are pulled-back to the background spacetime .Mo using a gauge choice X . In 
this sense, all variables treated below are tensor fields defined on the background 
spacetime Mq. We also denote the perturbative expansion of the pull-back of the 
variable Q on the physical spacetime M by 



Q 



e k \ k ' 



oo 

2 IfeijfeM 

k'k=0 



(k,k'\ 



'0- 



(3-7) 



as in Eq. (l2~231) . 

Once we have derived the formulae of the perturbative Riemann curvature (see 
flj.lj) of each order, it is straightforward to derive corresponding formulae of the Ricci 
curvature (* J3.2|) . scalar curvature (fc )3.3j) . Einstein tensor ( ffi.4|) . and Weyl curvature 
(i J3.5j) . We also derive the perturbative form of the divergence of an arbitrary tensor 
field of the second rank to check the perturbative Bianchi identities. 

3.1. Expansion of the inverse metric and the Riemann curvature 

Following the outline of the calculations explained above, we first calculate the 
perturbative expansion of the inverse metric. The expression for the inverse metric 
can be readily derived from the expansion (|2-31|) of the metric g a b and the definition 
of the inverse metric 

g ab g bc = S a c . (3-8) 

We also expand the inverse metric g ab in the form Q3-7JI . Then, each term of the 
expansion of the inverse metric is given by 



(l,l)-afe 



-^h a \ (2^) = (1,0), (0,1) (3-9) 

2 (|.f )h ac (f.f \ b - Mh a \ (p, q) = (2, 0), (0, 2) (3-10) 

(O.ljfcca (1,0)^6 + m h cb (1,0)^ a _ (341) 



To derive the formulae for the perturbative expansion of the Riemann curvature, 
we have to derive the formulae for the perturbative expansion of the tensor C c ab 
defined in Eq. (|3-5|) . The tensor C c ab is also expanded in the same form as Eq. (|3-7|) . 
The first-order perturbations of C c h have the well-known form 12 ' 



ab 



V7, (p,q)h , c 



-V C M hab 

2 



=: H. 



ab 



(3-12) 



where (p,q) = (1,0), (0,1), and in the brackets of the variable H ab c [( p ' 9 ^] 

indicates that H ab c [( p '^h\ is constituted of three covariant derivatives of the per- 
turbative metric ^ p ' q ^h ao . In terms of the tensor field defined by Eq. I|3-12|l . the 
second-order perturbations of C c ab are given by 



(p,q)rtc 
w ab 

(1,1W 
ab 



H ab' 



(1,%] _ ih0) h cd \{0,l) h ] _ {0,l) h ed Ul,0) h 



(3-13) 
, (3-14) 
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where (p, q) = (2, 0), (0, 2) in Eq. (l3~13|l . 

The Riemann curvature (|3-6j) on the physical spacetime A4 can be expanded in 
the form ()3-7j) . The forms of the perturbative Riemann curvature up to second order 
are given by 

lh abc 

(p,q)n , d 

±v abc 



-2Vr (P>l)(l d , 

-2V [o (p,9) C d 6]c + 4 (f-f )C e c[a (i -DC" 



(p,g) = (l,0),(0,l), (3-15) 



d b]e , (p,g) = (2,0), (0,2), (3-16) 

^R abc d = -2V [a (W qc + 2 (^)C e c[a (^)C rf 6]e + 2 ^C\ [a ^C%. (3-17) 

Substituting Eqs. (|3-12j) - ()3-14|) into Eqs. (|3-15)) - ([317)) . we obtain the perturbative 
form of the Riemann curvature in terms of the variables defined by Eq. (|3T2[) . This 
perturbative form of linear-order is simply given by the replacement 



(p,g)(7< 



fib 



H. 



ab 



(P,q) h 



(3-18) 



in Eq. ()3-15(l . On the other hand, the perturbative form of the 0(e 2 ) and 0(A 2 ) 
Riemann curvatures are derived from the perturbative form of O(eA). For these 
reasons, we only present the derivation of the perturbative form of O(eA) in terms 
of the variables defined by Eq. ()3-12j) , 



(1,1)5 



a be 



+2H, 



de 





H b]ce 


\o,i) h 


+ 2 V e 


\o,i) h 


H b]ce 





+2^%^V [a H b]ce 



(0,l) h + 2 (0,D^e V[a ^ 



(3-19) 



as the second-order perturbative curvature. 

To write down the perturbative curvatures (|3-15|) and (|3-19|) in terms of the 
gauge invariant and variant variables defined by Eqs. (|2-32j) - (|2-34j) . we first derive an 
expression for the tensor field H ab c in terms of the gauge invariant variables, 

and then, we derive a perturbative expression for the Riemann curvature. 

First, we consider the linear-order perturbation (|3-15[) of the Riemann curvature. 
Using the decomposition ()2-32(l and the identity R^ abc ^ d = 0, we can easily derive 
the relation 



Habc 



Habc 



where the variable H a bc p' 9 ^] is defined by 



H 



abc 



H ab C 



(p,«/)ft 

<p>*>H 



9cdH a 



1 



b) 



(3-20) 

(3-21) 
(3-22) 



Clearly, the variable H ab c [( p ' 9 ^] is gauge invariant. Taking the derivative of H a b c 
and using the Bianchi identity V[ a i? bc ] rfe = 0, we obtain 



(p,<?)e> 



abc 



- 2 ^[a H b]c ' 



+ £ 



(v,q)xH-abc i 



(3-23) 
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where (p,q) = (1,0), (0,1). 

Next, we consider the second-order curvature perturbation. We first consider the 
O(eA) metric perturbation as mentioned above. Inspecting the definition ()2-34j) of 
the gauge invariant variable of O(eX) metric perturbation, we first define the variable 



+ 2 {^O-MX^^X + £(0,l)x £( 1 >o)x) 9ab- 



(3-24) 



As in the case of linear order, we evaluate the tensor H 



ab 



and obtain 



2H ab C 



(i,i>£ 



2H. 



ab 



£(p,i)x \Hab° 

£{X,0) X [H ab C 

(o,i) h 



H a bd 



(o.i) h 

+ H a bd 



+ H ab ' 
+ H ab ' 



+ H r c 



ab 



(o,i)J 

£(o,i)x9 cd 
£(i,o)x9 



(o.D^e + (o,!)^ A (v a V 6 (1 ' 0) X d - i?^ 1 ' ^ 6 ) . (3-25) 



In the derivation of this expression, some formulae, which are summarized in Ap- 
pendix [XJ are useful. After straightforward calculations, we obtain 



^R abc d = -2V^, d \ {hl} H 



[a b]c 



+2H< 



de 



(1,0^ 



H,, 



+ 2H f 



de 



b]ce 

+ (0,l^d / £ {10)x R ahc e - ^->°>R abc e 



(1,0^ 



(1,0)5 d 1 /> p t 
n abc ~ 2 (lfi) X n abc 



(0,l) R 



1 



abc 



- £(o,i)x H-abc 



(m)# 



(3-26) 

, includes the gauge degree 



The first term on the right-hand side, — 2V r H a 

° 5 [a o\c 

of freedom, because the variable ^ l,1 ^H a b is transformed as a linear-order metric 
perturbation. Since we have already assume that the linear metric perturbation is 
decomposed as Eq. (|2-32j) . we can also decomposed the variable ^ l,l n~L a b as 



^>H ab =: ^>H ab + 2V {a ^X b) , 



4>), 



(3-27) 
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as pointed out in a previous paper. 8 ) The variables {1 > iy H a b and Wx b are gauge 
invariant and variant parts of O(eA) metric perturbation. Then, as in the case of 
linear order, we obtain 



2V r „ H, 



[a b]c 



(3-28) 



Further, using Eq. (|3-23|) . we reach the final form of the perturbative Riemann 
curvature of O(eA): 



abc 



-2V r # 



[a b]c 



+2H, 



de 



"(1,0)^" 


H b]ce 


"(0,1)^" 


+ 2 V e 


'(0% 


H b]ce 


'(1.0)^" 



+2( 1 '°)H e d V [a F % e + 2(°' 1 )H e d V [aJ Ff b]c 



+ £(0,l)x (1 '° ) - R a6c d + ^(i,0)x (0,1) ^, 



o 6c 



+ i?fi.i 



-£(i,o)x-£(o,l)x 



abc 



(3-29) 



The first three lines on the right-hand side of Eq. (|3-29|) are the gauge invariant part 
and the remaining two lines are the gauge variant part of the perturbative Riemann 
curvature of O(eA). 

The perturbative Riemann curvatures of 0(e 2 ) and 0(A 2 ) are simply obtained 
through the replacement of the variables in Eq. (|3-29f) of O(eA). To obtain the 
Riemann curvature of 0(e 2 ), we consider the replacements of the variables 



(0,l) X a ^ (l,0) x a (0,1^ _ (1,0)^ 



ab 



(3-30) 



in Eq. (|3-29|) . Similarly, to obtain the 0(A 2 ) Riemann curvature, we consider the 
replacements of the variables 



(l,0) x a ^ (0,l) x a j (1,0^ ^ (0,1)^ 



ab 



(3-31) 



in Eq. (|3-29j) . These replacements are consistent with the definitions (|2-33j) and 
(|2-34jl of the gauge invariant variables of 0(e 2 ) and 0(A 2 ). Hence, we obtain the 
perturbative forms of the Riemann curvatures of 0(e 2 ) and 0(A 2 ) as 



ip,l)v> d 

AL abc 



-2V r 



[a b]c 



+4H, 



de 



H. 



b]ce 



7? d 

^abc i 



(3-32) 



where (p,g) = (2,0), (0,2). 

Equations ()3-23|) . ()3-29(l and (|3-32|) show that all variables defined by 



"abc 



(p,q) R 



abc 



£ 



(p,i)X^abc 



(3-33) 
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for (p,q) = (1,0), (0,1), 
for (p,q) = (2,0), (0,2), and 

d n (1,0)S d P (0,1)S d 



£(p.i) X ~ £(P,P) X ) R a bc (3-34) 



afec 

A 1 .!)* ~~ o^fo.^X^t 1 . ^ ~~ 7;£( 1 ,°)x£( <Vx) Rabc (3-35) 



are gauge invariant. These indeed do have the same forms as the definitions ()2-35[) - 
(|2-37j) for each order gauge invariant variable for an arbitrary field. 

Here, we derive the perturbative formulae of the Riemann curvature R a bcd, which 
are used in the derivation of the Weyl curvature C abcd . For this purpose, we first 
expand the definition 

Rabcd = 9edRa.bc * • (3-36) 

The form of the perturbation of R a bcd at each order is derived from the formulae 

{P ' q) Rabcd = {P ' q) 9edR a bc & + ded^ Kbc \ (P, l) = (1,0), (0, 1), (3-37) 
( 1 ' 1 ^ = ( 1 ' 1 kA c e + (1 '°W ( °' 1) ^ c e 

+^gj i m abc e + gJ X ^R abc e ■ (3-38) 

The formulae for M R abcd with (p,q) = (2,0), (0,2) are derived using the replace- 
ments (|3-30j) and (|3-31|) of the perturbative variables as in the case of the Riemann 
curvature R abc d - The explicit form of each order ^ p,q ^R ab cd is as follows: 



M Rabcd = -2V [a H b]cd 
for (p,q) = (1,0), (0,1), 
{p ' 9) Rabcd = ~2V [a H b]cd 



+ (P<«>H d e R abce + £ Mx R abcd (3-39) 



b«>H + R abce ^>H d e + 4H d[a e b«>H 



H, 



b]ce 



abed 



(3-40) 



for (p,q) = (2,0), (0,2), and 



^R 



abed 



-2X7[ a H b ] cd 



{1 M + RabJ 1 ' 1 ^ 



+2H 



d[a 



H 



b\ce 



+ 2H, 



d[a 



H 



b\ce 



+£(i,o)x ^Rabcd + £(o,i)ir 



Wx 



labed 



(3-41) 



The perturbative forms (|3-39|) - ()3-41|) also show that the variables defined by 

iP ' q} Kabed ■= iP ' g) Rabed ~ £( M ) X Rabcd, (p, <?) = (1, 0), (0, 1), (342) 
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(P ' q} Ra b cd ■= M Rabcd ~ 2£ { P^ )x ^ P 2)R abcd 

- [£ Mx - £ 2 iM)x ) Rabcd, (P, q) = (2, 0), (0, 2), (3-43) 

^^R-abcd '■= ^'^Rabcd ~ £(o,i) x ^'^Rabcd ~ £ (i^x^'^ Rabcd 

£(i>Vx ~ 2^ {0 - 1) x£( 1 '°)x ~ 2^ ( - 1 ' 0) x£(°' 1 )x^j Rabcd (3-44) 
are gauge invariant. 

3.2. Ricci curvature 

Contracting the indices b and d in Eqs. (|3-23|) . ()3-29|) and ()3-32|) of the pertur- 
bative Riemann curvature, we can derive the formulae for the expansion of the Ricci 
curvature, 



Mr 



ab 



+ £ (p,q)xRab 



(3-45) 



for first order and 

M Rab : 



+ 4iJ, 



cd 



a 
E E\ 



H 



c]bd 



+^Wv [a H b] / (5.§)W 



MR, 



ab 



+2£ { Pv )x ^U)R ah +[£ 



(p,q)X ~ ) x ) R ab, 



(3-46) 



-^ [a H c]b 



+2H, 



cd 



'(1,0^' 


H c]bd 


"(0,l)ft 


+ 2^ 




H c]bd 





+2^Wy [a H c]h d [^>H] + 2(°'%/V [aJ ff c]6 d [WW 

+ <£(0,i)x ^ 1,0 ^o6 + "£(l.o)X ^ob 

C 1 ' 1 )* ~~ 2'^ (0 ' 1) ^'^ (1 ' 0)x ~ 2^ (1 ' 0)x ^ (0 ' 1)x ) Rab ' (3-47) 

for second order, where (p, g) = (1, 0), (0, 1) in Eq. (|3-45(l and (p, q) = (2, 0), (0, 2) in 
Eq. (EH|. 

It is trivial from the derivation that Eqs. IJH-45JI . (|3-46|) and (jH-47|) show that the 
variables defined by 

{M) Tlab = iP > q) Rab ~ £ Mx Rab (348) 

for (p,q) = (1,0), (0,1), 

^">Jl ab = { P^R ab - 2£ ( Pv ) /2^R ab - ^£ Mx - £%v )x ) Rab (3-49) 
for (p,q) = (2,0), (0,2), and 

( 1 ' 1 'TZah = ^'R-ab ~ £ (O^x^'^Rab ~ £ \i,o)x^° '^Rab 

£( 1 > 1 )x ~ 2~£(°> 1 )x£(- 1 '°)x ~ 2~£( 1 >o)x£(°, 1 )x J Rab (3-50) 
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are gauge invariant. These also have the same forms as those in the definitions (j2-35jl - 
(|2-37[) of each order gauge invariant variable for the perturbation of an arbitrary field. 



3.3. Scalar curvature 

The scalar curvature on the physical spacetime A4 is given by 

R := g ab R ab . 



(3-51) 



To obtain the perturbative form of the scalar curvature, the left-hand side of Eq. (|3-51|) 
is expanded in the form (|3-7|) and the right-hand side is expanded by using the Leib- 
niz rule. Then, the perturbative formula for the scalar curvature at each order is 
derived from perturbative form of the inverse metric (|3-9f) — llf) and the Ricci cur- 
vature (|3-45|) - ()3-47j) . Straightforward calculations lead to the expansion of the scalar 
curvature. 

Using (j3-9jl and (|3-45j) . we obtain the first-order perturbative form of the scalar 
curvature as 



(p,q) R 



-2V [a ff 6] ' 



(™W| - R ab ^>H ab + £ Mx R, 



(3-52) 



where (p,q) = (0,1), (1,0). Then, using (32J), (ETTTll . (ET25jl and (|3~T7|) . the pertur- 
bative scalar curvature of 0(e\) is found to be given by 



ab 



'C 1 ' 1 ^] + R ab (2( 1 '°^ c(a (°' 1 % c - ^>H ab ) 





tt a 

H c\ d 




+ 2^ [a cd 




tt a 

H c] d 


"(0,l)ft 



+2H [a cd 

+ 2^°>H c b V [a H b] ac l^hi] + 2^ 1 >H c b V [a H b] ac \ (1 >°>H 
+2^°>H ab V [a H d]b d [^>H] +2(°' 1 )H a6 V [a F d]fe rf [^°>H 

+£ (0.1)V 

-£(o,i) X £(i,o) X ■ .-, 



+ ( £(i,i)x ~ oAo,i)X"£(i,o)y - t;£( 1 <o)x£(o,i)x I 



(3-53) 



To derive the perturbative scalar curvature of 0(e 2 ), the replacement ()3-30j) of 
the variables is applied to Eq. (|3-53|) , Similarly, the replacement ()3-31j) is applied 
to Eq. (|3-53j) when we derive the perturbative scalar curvature of 0(X 2 ). Then we 
obtain the perturbative form of the scalar curvatures of 0(e 2 ) and 0(\ 2 ): 



[p,q) R 



ab 



cd 



a 



H 



c] d 



^>H\+m>H b V [a H b] ac m>H 



+4 P 2> P 2>H ab V [a H d]b d $4>H 



'Wx 

where (p,q) = (2,0), (0,2). 



+2£ ( p .P)„ 2,2 'R + ( £(p,q) X — °^(f,fw 



R, 



(3-54) 
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It is also trivial from the derivation that Eqs. (|3-52|) - l|3-54|) show that the vari- 
ables defined by 



(v^yji = Mr _ £ {Ptq)x R 

for (p,q) = (1,0), (0,1), 

(p,^ = Mr _ 2 £ {M)j flm - (£ Mx - £^ 

for (p,q) = (2,0), (0,2), and 

(i,i)^ = (i,Dr _ £{0A)x M)r - £ (U))X ^R 



R 



(3-55) 



(3-56) 



R 



(3-57) 



are gauge invariant. These too have the same forms as those given in the defini- 
tions ()2-35j) - (|2-37j) of each order gauge invariant variable for the perturbations of an 
arbitrary field. 

3.4. Einstein tensor 

Next, we consider the perturbative form of the Einstein tensor. The Einstein 
tensor on the physical spacetime M is defined by 



G, 



ab 



Rab — -^QabR- 



(3-58) 



As in the case of the scalar curvature, the left-hand side of Eq. (|3-58j) is expanded 
in the form ()3-7(l , and the second term on right-hand side of Eq. (J3-58|) is expanded 
by using the Leibniz rule. Then, the perturbative formula for the Einstein tensor at 
each order is derived from the perturbative form of the metric (j2-32|) - H2-34|) . those 
of the Ricci curvature l|3-45|) - (|3-47|) . and those of the Ricci scalar l)3-52Jt - (|3-54J) . 
The linear order Einstein tensor is given by 



ab 



- 2 ^[a n d\b 
~ ! <~-£(p,q)X^'abj 



+ 9ab^\ c K 



[c' L d] 



cd 



<P«>H - R (P^>H ab + -g ab R cd ^>H cd 



(3-59) 



where (p,q) = (1,0), (0,1). Next, using Eqs. (t2~^|) - (l2^4l . (l3~4T]l . (|3~52l and (l3~53l . 

the Einstein tensor of 0(e\) is found to be given by 



(1,1)6, 



ab 



-2V [a ff c]6 c 



(1,1)^ 



+2H, 



tic 



"(!,(%" 


Hd]be 


"(0,1)^" 


+ 2 V e 




Hd]be 


'(1,0^ 



+2( 1 '°Vv [o H d]6 e ^>H +2^>H e d V [a H d]h 



~g ab (-2V [c H d] cd 



^>H + R de (2^ x >H c d ^°>H ec - 



+2H, 



de 



(1,1% 



H 



d] e 



(O.Dft 



2H, 



de 



(0,1)^ 



H 



d] e 



Second-order Gauge Invariant Perturbation Theory 



21 



+2^>H e d V [c H d] ce W>H\ +2(°> 1 >ft e d V [c tf d] ce \^°>H 



[c d]e 



+ £(0,i) x ( * 1 '°* > Gab + £(i,0)x ( ' ' 1 * > Gab 



(3-60) 



Then, through the replacements (|3-30|) and (|3-31|) . the perturbative forms of the 
Einstein tensor of 0(e 2 ) and 0(A 2 ) are given by 



(p,i)q 



ab 



1 

~2 9a 



cd 



H. 



c]bd 



cd 



de 



H d] C e 



'E £ 



+2£ (f.§)x (f ' § '^ 



^ 2 ' 2 ^Gab + ( £( P ,q)x ~ £%,§) x 



G a b, 



(3-61) 



where (p,g) = (2,0), (0,2). 

Further, we also derive the formulae for the perturbation of the Einstein tensor 



G a ° := g bc G 



(3-62) 



on Ai. Because the derivation of these formulae is similar to the above perturbative 
curvatures, we only present the final results: 



P C k 

(p>q)X a 



(p,q) G< b = iXigb (p,^ 
for (p,q) = (0,1), (1,0), 

(p,q)Q a b = {\)gb [(P,«) W 1 + i?)gb \(P,Q)f{ } 



(3-63) 



(3-64) 



22 



K. Nakamura 



for (p,q) = (0,2), (2,0), and 

(1,1)^6 = {l)gb r(l,l) W l + (2)^6 |'(l,0) Wj (0, 1 ) W 



where 



+^ a b (2V [eJ ff dl ed [yl] + J R ed ^ 



-2V [a i^ bd [A]-A cb R a 

( 2 V [A, 5] := 2i? a ,B c (^) c + 2^/ e [A] [-B] + 2flj tt * [5] i/ d] 6 e [A] 

+2A e d V [a ^ fee [5] + 2B/V [a H d] 6e [A] 
+2A c V [a ^ cd [5] + 2£? c 6 V [a tf d] cd [A] 



(3-66) 



-S a ° [R de B f {d A^f + H [f «* [A] H/ e [B] + H [f «* [B] H^ e [A] 

+2A e d V [f H d] [fe] [B] + 2B e d V [f H d] [/e] [A]) . (3-67) 

We note that ^G a b [*] and ^Q a b [*, *] in Eqs. (|3-63|) - (|3-65|) are the gauge invariant 
parts of the perturbative Einstein tensors, and each expression (|3-63|) - ()3-65|) has a 
form similar to Eqs. (|2-35j) - ()2-37j) . respectively. 

3.5. Weyl curvature 

Here, we consider a perturbation of the Weyl curvature, which is useful to study 
some physical situations. In m-dimensional spacetime, the Weyl curvature is defined 
by 



Cabcd '■— R, 



abed 



m — 2 



{9a[c R d]b - 9b[c R d}a) + 



(m — l)(m 



-R9a[e9d]b- (3-68) 



Using the perturbative formulae for each order perturbation of the Riemann curva- 
ture (pH^ - iPRTj) . of the Ricci curvature (|3~15)l - (i3^oTI)l . of scalar curvature (j3~5oT) 
(|3-57[) . and of the metric perturbation H2-32j) - (|2-34[) . we can explicitly confirm that 
the perturbative forms of the Weyl curvature at each order are given by 



^P^Cabcd — ^'^Cabcd + £(p,q)xCabcd 

ip ' q) C abcd = M C, 



for (p,g) = (0,l),(l,0), (3-69) 



abed + 2£ ( P,9) x ^ 2 ' 2 ^C afecrf + ( £( 



f 2 



C, 



abed 



for (p,q) = (0,2), (2,0), (3-70) 

^'^Cabcd = ^'^Cabcd + £(0,i)X ( ' 1 ' ^Cabcd + £(h°)X ^'^Gabcd 

+ (£(i,i)x ~ ^(i&X'Z^x ~ 2^(°' 1 )-x:^'C 1 ' )x^ C a bcd, (3-71) 
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where 



M C ab cd = {P ' q} Ra b c d ~ { (P ' q) Ha[cR d ] b ~ {P ' q} H b lcRd]a 



+9a[c {P ' q} R- d]b -9 b[ c {M) n d] a} 

(rn-l)(rn-2) { iP ' q}R9 ^ + R 



+R g 



'a[c iP ' q) H d ] b } 



(3-72) 



for (p,q) = (0,1), (1,0), 

(P'l^Cabcd = ^ P ' q ^T^ahcd ~ 



m-2 



"(m- l)(m-2) 



■abed 



Xp ' q} Ha[cR d] b - {p ' q} Hb[cR d ] a 

+9a[c iP ' q} K d ]b-g b[ c iP ' q) K d ]a 

{P '^5a[ C ^] fe + « (P ' 9) Wa[ C ^]6 + Rg a[ c iP ' q} H d ]b 

+2^>Rg a[c ^'^b + 2^i)R ^>H a[cgd]b 
+2R ^i>H a[c { ^>H d]b \ (3-73) 

for (p,q) = (0,2), (2,0), and 

(1,1) ?4[ci?d]& - (1,1) H6[ C -Rd]a 

-^^[c (1 '%a " (1 ' 0) H6[c (0>1 V 
(1,1) K<?a[c<7# + R (1 ' 1} Ha[c9 d ]b + R9a[c (1 ' 1} Hd]b 

+^>JZ 9a[c ^>H d]b + ^>R^>H a[cgd]b 
+^°>R^ki a[c g d]b + ^>JZg a[c ^>H d]b 
+R ^°>H a[c ^>H d]b + R ^>H a[c ^>H d]b } .(3-74) 

We also derive the formulae for the perturbative expansion of C abc d from 



^Cabc d = 



m-2 



(m - l)(m - 2) 



^abc — 9 ^abce 



(3-75) 



Since the Weyl curvature is traceless, i.e., C ' , d = 0, we can also verify the formulae 
derived here by confirming this traceless property. 

Actually, we explicitly confirm this traceless property with the following formulae 
for the perturbative Weyl tensor at each order: 

M C ab / = M C a bc d + £(^xC abc d for (p, «/) = (0,1), (1,0), (3-76) 
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where 



M C abc d = M C abc d + ,i •* } C abc d + (£ Mx - £ 2 iU)x ) C abc d 

for (p,g) = (0,2),(2,0), (3-77) 

C abc - L abc + ^(O.l)X L abc + ■*< W)X C abc 

+ y£a,i)x ~ 2 ■Ai,o)xAo,i)x ~~ 2^ (0 ' 1) x^ (1,0) -f j Cafec^' (3-78) 

(P ' 9) Ca bc d == - iP ' q} H de C abce + ^^Cabce, (3-79) 

_2(f-§>H de ( i'i)C afcce , (3-80) 
(i-i) Ca6c d = - (^>H de - ^°>H f e - 0"°>H f d (0 '% e/ ) C abce 

+9 de ^Cabcc ~ ilfi) Cabcc - ^C abce . (3-81) 

By using the fact that C abc b = for the Weyl curvature on the background A4o, 
Eqs. (|3~72l . (|3~4"5i (|3~521) and ((2123), a straightforward calculation yields 

{p ' q) C ab c = (3-82) 

for (p, q) = (0, 1), (1, 0). Further, with straightforward calculations using Eqs. (|3-74|) . 
(ETT21) . (f3~T77|l . (jHSSJ), dHUDJ), (E3S1), (iHTTTl and (HUH), we can explicitly confirm the 
identity 

= 0, (3-83) 

and hence, we obtain 

{p ' q Kbc = (3-84) 

for (p, q) = (0, 2), (2, 0) through the replacements (|3-3()|) and (|3-31|) of the perturba- 
tive variables. The property ^ p ' qS> C abc b = is trivial from the definition of the gauge 
invariant part of the Weyl curvature. However, this trivial result gives us great 
confidence in the formulae derived here. 

3.6. Divergence of an arbitrary tensor of second rank and the Bianchi identity 

Here, we consider the perturbation of the Bianchi identity and the divergence 
of the energy momentum tensor, which are derived from the divergence V a T fe a of 
an arbitrary tensor field T b a of second rank. The operator V a are the covariant 
derivative associated with the metric g ab on the physical spacetime M. . As discussed 
above, V a is pulled back to the background spacetime Aio as the derivative operator 
X*V a (X- 1 )* by choosing a gauge X. Further, the operation of V a as an operator 
on Mq is represented by the covariant derivative V a , which is associated with the 
background metric g ab on .Mo, and the tensor field C c ab defined by Eq. (|3-5|) , Hence, 
we may concentrate on the Taylor expansion of the equation 

V a f b a = V a f b a + C a ca f b c - C\ a f c a (3-85) 
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to derive the perturbative form of the divergence of an arbitrary tensor field of second 
rank. The tensor field T a b , which is pulled back from the physical spacetime Ai € to 
the background spacetime Mq, is expanded as Eq. ()3-7|) , Following the definitions 
(|2-35j) - (|2-37j) of gauge invariant variables, the gauge invariant variables (P' 9 ^ a for 
each order perturbation (P ,q>r T a b are defined by 



{p,q}j a ._ {p,q)rr a _ n rp a 



(p,9) = (l,0),(0,l), (3-86) 



c*.4>x 

-{£(v, q)x -£ 2 { v, i)x }T b a , (p,g) = (2,0),(0,2), (3-87) 

v ; T 6 := ^ >L b - £(i,o)x K r h ~ *(o,i)x a b 



£(^)x ~ ^^^x^^^x ~ 9 ■Zw-Vx^wix } T b a . 



We also expand 



k',k=0 

3.6.1. Linear order 

A simple expansion of Eq. (|3-85j) yields 

ip ' q \v a f b a ) = V a ^ q) f b a + (H ca a \&<i)h\ + V c V a Mx a ) T t 



(3-88) 



(3-89) 



H, 



bac 



(P' q >H\ + V 6 V a b>9>X c + R acb e ^X e ) T ca (3-90) 



for (p, q) = (0, 1), (1, 0). Then, using the gauge invariant variable defined by Eq. 
we obtain 



^ q \v a f b a ) = V a ( p ' qy T b a + H ca a [( p ' q >H 

+ £ d>,q)x^ 'a,T b ° 



T, 



H, 



ba 



<p>*>H 



rp a 



(3-91) 



for linear order [(p, q) = (1, 0), (0, 1)], where we have used the formula (|A-4[) . 

Now, we check the linear-order perturbation of the Bianchi identity by using 
Eq. ()3-91j) . To do this, we regard T a b as the Einstein tensor G a b . Further, the gauge 
invariant part of the tensor T a b is regarded as the gauge invariant part of the linear 
order Einstein tensor: 



rp b s~1 b 

(p.9>r b _» 0)n b \(P'i}]-/ 
Through these replacements, we obtain 



(3-92) 
(3-93) 



M [v a G b a ) = v a ^g a b [^ q >H 

J r£(v,ci)x'^aG b a . 



+ H ca a 



& q >H 



Gi, 



H ba C 



G c a 



(3-94) 
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On the other hand, a direct calculation using Eq. (j3-66|) yields 



v a (1) g h a \ {p ' q) n 



Gh° + H ba c 



(3-95) 



as an identity. Therefore, the linear-order expansion of the divergence of the Einstein 
tensor is given by 

(M) (V a (V) = £i P , )x V a G b a , (p,q) = (1,0), (0, 1). (3-96) 

Because X7 a G b a = for an arbitrary spacetime, we can easily see that {p ' q \V a G b a ) = 
0, identically, at linear order. 

3.6.2. Second order 

Next, we consider the O(eA) perturbation of V a T b a - A straightforward calcula- 
tions yields 

H C ad 



+ H 



l bad 



ba 



-H 



ba 



(°>% ^>H da + H cad &°>H ^>H da -H ca a ^>H 
(0,1^ (1,0^0 + Hbad J(i, 0)H j (o.Dftdc _ Hba c ^(1,1)^ 

(l.O^j {0,iy Tc a + Hca a Jd.O^j (0,1)^ c 
(0,1)J (1,0% a + fl a (0,1)^ (1,0% c 



T,, 



J c 



(V„V) 



(3-97) 



Applying the replacements ()3-30j) and (|3-31|) of the perturbative variables to Eq. ()3-97(l , 

the 0(e 2 ) and 0(A 2 ) perturbations of the divergence of a tensor T a b are obtained as 



— {%H ca d 
+ ( ^H ba d 



m a 



-2F, 



ba 



(M>h\ ^K a + 2H C /[U>H\ (M>T b c 
+2£ (U)x (U\v a f b a ) + {£ Mx - 4f,*) x } (V«V) ,(3-98) 

where (p,g) = (2,0), (0,2). 

As in the linear-order case, we can also check the Bianchi identities of the second 
order perturbations using Eqs. (|3-98|) and ()3-97j) . To do this, we regard T a b as the 
Einstein tensor G a b and apply the following replacement: 



r 1 b 
Ll a > 

(% b 
=>a 



for (p,g) = (l,0),(0,l), 



(3-99) 
(3-100) 
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for (p,g) = (2,0), (0,2), 
First, we consider the perturbative Bianchi identity of O(eA): 

^\v a G b a ) = v a («s 6 ° p 1 ^] +( 2 )g 6 a [c 1 ^,^ 1 ^]) 

- (H cad [< >%] &°>H* + F cad [(^)w] - tf ca a 

(i.o^j (o^a + ^-po^j (0,1)^ c 

+£ (1)0)x ^'^(VaG^) + £ C o,Dx (1 ' 0) (V a G b a ) 
Here, we note that the identity (|3-95|) implies 



(3-101) 
(3-102) 



G, 



G c a 



~ H ba C 



-H ba c 



- ~ H ca ° 



"(1,1)^ 



G fe c + if, 



6a 



G c a . 



(3-103) 



(3-104) 



Further, straightforward calculations lead to the following identity: 



tt a 

~ H ca 


"(1,0)^" 




'(O.Dft" 


i_r a 




(1 V 


(1,0)^" 


+ H ba C 


"(1,0)^" 


(Dg e a 


"(0,l)ft 


+ ^ba £ 






"(1.0^" 



(Hbad 


(0,1)^" 


(l.O^dc 




"(!,<%" 








(H ca d 




(1,0W 




"(1,0^ 


(0,l>^arf N 


)G b c . 


(3-105) 



Using the identities Q3-1U4J1 and (|3-105|) , we easily find 

^\V a G b a )= £ (mx ^\v a G b a ) + £ i0 , 1)x W>(y a Gf) 

}iy a G b a ). (3-106) 

Since the Bianchi identities V a G b a = on the background spacetime A4q and those 
of the linear-order perturbations (j> ' q \V a G b a ) = [for (p,q) = (1,0), (0,1)] have 
already been confirmed, we have also confirmed the identity 

(^(VaG^) = 0. (3-107) 



28 



K. Nakamura 



Applying the replacements 1)3- 30 j) and (|3-31|) of the perturbative variables to 
Eq. (|3-103|) . we obtain the 0(e 2 ) and 0(A 2 ) Bianchi identities: 



-( 



2H, 



cad 



+ 2H bl 



■ad 



-2iT, 



6a 



+ 2H r 



G b c 
G c a 



(£,«) 

= 0, (3-108) 
where (p,g) = (0,2), (2,0). 

3.7. Einstein equations 

Finally, we consider perturbations of the Einstein equation at each order. First, 
we expand the energy momentum tensor as Eq. (|3-7jl and impose the perturbed 
Einstein equation of each order, 

(p,o)Qb = 87T q (p,q) Ta b ( 3 . 109 ) 

Then, we define the gauge invariant variable ^ p,q ^T a b for the perturbative energy mo- 
mentum tensor at each order by Eqs. (|3-86[) - (|3-88j) . Then, the perturbative Einstein 
equation at each order is given by 



for linear order [(p, q) = (0, 1), (1, 0)] and 

(Dg o & r (1) l) H l + (2)g a b [(1,0)^(0,1)^1 =87rG (hl) Ta \ 



^>a 



ip,D n 



+ {2) Qa h 



(U) n (f.f)w =8ttG <**)t 



(3-110) 

(3-111) 
(3-112) 



for second order [(p,q) = (0, 2), (2, 0)]. These explicitly show that, order by order, 
the Einstein equations are necessarily obtained in terms of gauge invariant variables 
only. 



§4. Summary and Discussions 



In summary, we have derived some formulae for second-order gauge invariant 
perturbations, namely those of the Riemann, Ricci, scalar, Einstein, and Weyl cur- 
vature tensors. We also derived the formulae for the divergence of an arbitrary tensor 
field of second rank for perturbations at each order. These perturbative curvatures 
have the same forms as those in the definitions (|2-35|) - (|2-37|) of the gauge invariant 
variables for arbitrary perturbative fields which are proposed in a previous paper. 8 ) 
These are useful for investigating physical problems using perturbation theory. 
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Through the derivation of these formulae, we have confirmed some facts. First, 
if linear-order gauge invariant perturbation theory is well established, its extension 
to higher orders and multi-parameter perturbations is straightforward. Second, the 
perturbative Weyl curvature at each order preserve the property of the tracelessness 
of the Weyl curvature. Third, the perturbative Einstein curvature of each order 
satisfy the perturbative Bianchi identities at each order. Fourth, the perturbative 
Einstein equations at each order necessarily take in gauge invariant forms. These 
properties we have confirmed are trivial from their derivations. In particular, the 
fourth result is trivial, because any equation can be written in a form in which the 
right-hand side is equal to "zero" in any gauge. This "zero" is gauge invariant. 
Hence, the left-hand side should be gauge invariant. However, we emphasize that 
these trivial results imply that the formulae derived here and the framework of 
higher-order perturbation theory applied here are mathematically consistent at this 
level. 

Further, we note that in our framework, we specify nothing about the back- 
ground spacetime nor about the physical meaning of the parameters for the per- 
turbations. Our framework is based only on general covariance. For this reason, 
this framework is applicable to any theory in which general covariance is imposed, 
and thus it has very many applications. Actually, we are planning to apply it to 
some physical problems. The following are candidates of the physical situations to 
which the second order perturbation theory should be applied: the radiation reac- 
tion problem in gravitational wave emission; 13 ) stationary axisymmetric ideal MHD 
flow around a black hole or a star; 14 ) the correspondence between observables in 
experiment and gauge invariant variables; dynamics of gravitating membranes (for 
example, topological defects, 4 ) brane world, 15 ) and so on); perturbations of a com- 
pact star with rotation and pulsation; 6 ) Post-Minkowski expansion alternative to 
post-Newtonian expansion; 16 ) higher-order cosmological perturbations and primor- 
dial non-Gaussianity. 17 ) 

In particular, the gauge invariant form of the second-order perturbation of the 
divergence of the energy momentum tensor should be useful in considering the gauge 
problem in the radiation reaction of gravitational wave emission. In astrophysical 
contexts, it is natural to consider the situation in which a solar mass object falls into 
a supermassive black hole of mass ~ 10 6 M Q . This is one of the target phenomena of 
the observation of gravitational wave by LISA (laser interferometer space antenna for 
gravitational wave measurements). 18 ) In such a situation, the perturbation parame- 
ter is the ratio of the mass of the compact object to that of the central supermassive 
black hole. We may regard the second order perturbations as describing the radia- 
tion reaction effect of gravitational wave emission. By applying the gauge invariant 
formulation discussed here, we can exclude gauge freedom completely and thus there 
is no gauge ambiguities in the results. It would be quite interesting to apply the for- 
mulation discussed here to this radiation reaction problem. We leave this application 
as a future work. 

In addition to the radiation reaction problem of gravitational wave emission, 
there are many applications, some of which are listed above. Indeed, because there 
are few assumptions in our treatment, it is natural to expect that there are many 
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applications that are not listed above. We also expect that the formulae derived here 
will be found to be very powerful tools in many applications. 
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Appendix A 

Useful formulae 



In the calculations in the main text, the knowledge of the commutation relations 
of the covariant derivative and Lie derivative are useful. Here, we summarize the 
commutation relations which are used in the derivation in the main text: 



V a £xh = £xV a t b + X c R acb % + t c V a V b X c , (A-l) 

e± , yd t> e 
adb l ec "I" -A U adc 



^ a £xhc — £xV a tbc + X d R adb e t ec + X d R adc e tb. 



+t dc V a V b X d + t bd V a V c X d , (A-2) 

^a£xh C = £x^ ah" + XdR adb et e C - XdR ade% e 

+t/V a V b X d - t b d V a V d X c , (A-3) 
V a £xt b cd = £xV a tbcd + X e R aeb ■ f tf cd + X e R aec ft b f d + X e R ae J t bc f 

+t ecd V a V b X e + t bed V a V c X e + t bce V a V d X e , (A-4) 
V a £ x t bc d = £xV a t bc d + X e R aeb f t fc d + X e R a Jt bf d - X e R aef \J 

+t ec d V a V b X e + t be d V a V c X e - t bc e V a V e X d . (A-5) 
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